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Abstract 

The niinimal dark matter model is given a supersymmetric extension. A super SU{2)l quintuplet 
(— I ' is introduced with its fermionic neutral component still being the dark matter, the dark matter 

particle mass is about 19.7 GeV. Mass splitting among the quintplet due to supersymmetry particles 



is found to be negligibly small compared to the electroweak corrections. Other properties of this 
supersymmetry model are studied, it has the solutions to the PAMELA and Fermi-LAT anomaly. 
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^ ■ the predictions in higher energies need further experimental data to verify. 
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I. INTRODUCTION 

It has been known from many astrophysical measurements that the universe contains 
enormous of invisible, non-baryonic dark matter (DM) which is not included in the Standard 
Model (SM). Among various hypotheses for the nature of the DM, that of weakly interacting 
massive particles (WIMPs) is very attractive. Focussing on the DM problem, one can explore 
a simple WIMP model, that is the minimal dark matter model (MDM) [l|, |2|: adding to 
the SM a single matter X without introducing any additional discrete symmetry, and X is 
in a high dimensional representation of the usual SM SU{2)l x U{1)y electroweak (EW) 
symmetry. The stability of the DM candidate is guaranteed by the SM gauge symmetry and 
by the renormalizability. The minimality of the model lies in the fact that the new physics 
is determined by only one parameter, namely the mass M of the X multiplet. Therefore, 
the MDM is remarkably predictive. There are some extensions to the MDM |3|. 

As far as the particle physics is concerned, the SM provides a successfully description of 
presently known phenomena. However it will have to be extended to describe physics at 
higher energies. That often results in the gauge hierarchy problem. Supersymmetry (SUSY) 



u. 



offers a solution to this problem 

In this work, we make SUSY extension to the MDM. Note that in the so-called minimal 
SUSY extension of the SM (MSSM), the DM candidate, that is the lightest SUSY particle, 
is there only after introducing an extra discrete symmetry by hand, which is the i?-parity. 
In our SUSY MDM (SMDM), instead, we still follow the logic of MDM, the existence of the 
DM lies in the fact that the DM is in a high dimensional representation of the SM gauge 
group without using discrete symmetries. 

In section II, the SMDM is constructed. In section III, mass splitting of the X multiplet, 
the DM relic density, direct and indirect detection signatures of the SMDM are calculated. 
In section IV, the conclusion is made. In the Appendix, we give basic facts about the 
representation of the SU(2) group. 

II. SMDM 

The SMDM is simply constructed by supersymmetric extension to the MDM. The particle 
content is, in addition to that of the MSSM, the fermionic SU{2)l 5-plet X of the MDM and 



(1) 



its superpartner which a complex scalar 5-plet X. The charged components X'^ are slightly 
heavier than the neutral one X" due to quantum corrections [l|, y] ; and the superpartner X 
because of its soft mass Mgoft, is also heavier than X. Both X*^ and X will decay into X*^. 
The relic particle in the SMDM is still X° as in the MDM. The new parameters are M and 
Msofti the model is still predictive. As for the Lagrangian, in addition to that of the MSSM, 
we have 

C-SMDM = ^(Xt^a'^S^X, + Xt'a^'S^X,) - '-g^A'^.iX^'a^iT^^X^ - Xt,a^(T")/X,) 

\/2 ~ ^ ^ _ _ ^ 

- — ^2(X"(T")i^XjA'' + At"XtX^'')i^Xi - X"(T'^)/X^A" - At"X^*(T")/Xj) 

+ 1(Z}^X*Z}^X - M^ixn - ^M(X,X, + Xt^Xt^ + ]^g^D^X\T'^)i% 

- ^M^ „|XP. 

2 sojt\ I 

The component field notation has been used. In eq.(l), T"'s are generators of the SU{2)l 
in n representation. Xj and X* consist of the left-hand pairs of Xj, transforming in SU{2)l 
5 representation with the generator {T"')i^ and the complex conjugate representation with 
the generator iT"")*^ = (T'^)^*, respectively. They are not independent. Actually they are 
dual to each other under the SU{2)l- We write the Lagrangian in the form of eq.(Il]) just 
for convenience. Their superpartners compose a bosonic SU{2)l 5-plets X*. Both Xj and 
X* have trival SU{3)c x U{1)y quantum numbers (1,0). 

In terms of 4-component notation, we can define following spinors, 

x|/+2= I ^^' I ,^"'= I ^ ' I ={^^Y, (2) 




^+ = - .J,^" = - .J =(^+T, (3) 



(4) 



The neutral field \1'° is a Majorana field. 
The superpotential takes the simple from: 



W = Wmssm + ^MX^ (5) 



which gives X and X the same unbroken supersymmetry mass M . X get a soft mass after 
supersymmetry breaking. The D-term contribution to the scalar potential is: 



Vd = \9l 



1 '^ 

^ 0*t"0 + -X*T"X 



(6) 
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where denotes the SU{2)l scalars in the MSSM. Compared with the MSSM, the extra 
term is: 

(7) 
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2^2 



^(0*t»(X*T"X) + ( -X*T"X 



These couplings do not cause X dacay but annihilation into MSSM SU{2)l scalars, which 
give an extra negligible mass splitting between X*. 

Considering non-renormalizable terms of the lagrangian, there are dimension 5 operators 
X*-^'^'0j0j0fc0i/A for the complex scalar 5-plet X, and dimension 6 operators X*-''^'-?/'j0j0fc0;/A^ 
for the fermionic 5-plet X allowed by the SU{2)l x f/(l)y gauge symmetry, where ipi is the 
left-hand leptons or the higgsinos in the MSSM, eg. XH^HdHuH^/A, XHuHaHdH^/A, 
XH^HdHdH*JK\ XLHuHdHjA^ ,etc. 

We can generate these couplings by adding the corresponding higher dimension superpo- 
tential, eg. 

'^ non—ren 7"^ ~r 7"^ ~r ... , yo) 

the equation of motion for the auxiliary fields are: 

Fh, = -(^) = -il^H. + ^. + ^^^ + •••) , ^^^ 

This generates dim6 couplings for the fermionic 5-plet : XHuHdHaH^ / A^ , XLHuHdHu/ A^ 
where A k. lO^^GeV. These operators can induce 4-body decays with a typical life-time 
r ~ A^TeV~^ ~ lO^^s which is longer than the age of the universe (~ lO^'^'^s). So these 
couplings have no influence on the observed stability of the DM candidates. The F-term 
also generates the scalar dim. 5 operators but which are not suppressed by one power of A but 
two powers: fiXHuHdHdH*^ / A'^ , fiXLHdHdH^/ A"^ . The typical life-time of these operators is: 
T ~ A^TeV^^/i~^ ~ lO^^s which is long enough. Of course there are even higher dimension 
couplings of the complex scalar X, eg. XX*HuH*HdH*HdH*/A\ XX*LHuHdH*H*H*/A\ 
etc. These higher dimension operators can be neglected in considering the decay. 



Therefore the new introduced particles (X^^jX^^jX*^) only decay into (X^^,X^^,X°) via 
gauuge interactions. (X=''^,X=''^, X°) are quite stable. We will further study mass splitting 
among them in the next section, and see that the DM candidate is still X°. 

III. PROPERTIES OF THE SMDM 

The DM candidate in the SMDM model is still X°. Mass splitting due to SUSY particles 
is small because of SUSY breaking. 

A. Mass splitting 

Mass splitting should be studied in detail, like that in the MDM, because it can be 

calculated with little uncertainties in this simple model. Because of EW symmetry breaking, 

;he gauge kinetic terms gives the fermonic 5-plet X a mass splitting through loop crrections 

l|, bl, AM^y^r = iVf^ - M° ?^ Q2 X iqq MeV, where M'^ and M° are the pole masses of X'? 

and X°, respectively. 

The scalar particles X's also contribute to mass splitting of X. They are heavier than 
the fermions X's by a soft mass Mgoft which is generally expected to be about 100 GeV-1 
TeV. This further mass splitting is calculated by using the supersymmetric kinetic term, the 
second line of the eq.(l), at the loop level which involves X's and the gauginos. By using 
the two-component notation for fermions, the one-loop pole mass is written as (8| 

M%sY = Mil + ^Si« + 1sh«) (10) 

where Si'^,S/j'^ are the IPI self-energy functions as shown in Fig.l in Q = 0, +1 cases. 
In the diagrams of Fig.l, we denote the correspondent superpartners : 

0+2 ^ x+2, 0'-2 = X+2, 0+1 = X+\ 0'"i = X+\ 0° = X+\ (11) 

w+, ws are the superpartners of the SU{2)l gauge bosons and w^ = Vn^C^ + Vi2^C2, 
w~ = U^iCi + U12C2, W3 = N2i^Ni, i=l-4. (7+, C^ and iVj are the charginos and 
neutralinos of the MSSM. U, V, N are the unitary matrices diagonalizing the mass matrices 
of charginos and neutralinos p] . 

It is worthy to note that the T^d'^ term which is related to the B-term B^^XiXj, may also 
appear in the pole mass formula, and it does not cause divergences. But in our calculation 
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FIG. 1: One-loop corrections to the IPI self-energy functions to the Q = 0, +1 components of the 
SMDM. 

we do not consider it for simplicity, it is enough for us to break the supersymmetry only 
through the soft mass Mgoft- 

Using the superpartner notation mentioned in Sect. I we get : 



^t = T^[v;*iV^iifii(c'i,0+') + ^*iV^2i5i(c'2,0+') + -{u;,UiiBi{Ci, 



IGtt 



1 



+ f/2*if/2l5i(C2,0°)) + -(iV*2iVi25i(iVi,0+l) + iV*2iV225i(iV2,0+l) 
+ N;^Ns2Bi{Ns, 0+1) + iV*2iV425i(iV4, 0+1)], 



(12) 



SV = T^[f/rif^ii5i(C'i,0'-') + U;,U2lB^iC2A'~') + ^iV,\V^,B^{C,, -^^^ 



-•R 



167r2 



3 

2' 



1 



+ V,\V2lB^iC2, 0°)) + -(iV*2iVi25i(iVi, -1) + iV*2iV225i(iV2, 



'-U 



(13) 



+ f/*if/2l5l(C2,0-l))], 



(14) 



(15) 



+ \/;i\/2i5i(C'2,0'+^))], 

where Bi is the one rank two point integral 

Bi{p ,1711,1712) = -—^ —5 (16) 



]_ Ao(mi) - Ao(m2) + {ml - m l - p^)Bq{p^, mi, ms) 
2e 2p2 



with Aq and i?o being the Passarino-Veltman functions. 

Because all the superpartners (jf have the same mass M + Mgoft, we can simplify the 
above four equations to get the final result of the mass splitting due to SUSY particles : 

AM'3 - ^rvO _l vO _ vO _ y^O ^ 

= T7r^[-(^n^ii + U^lUll)B^iC^,(l>) - (^l^si + U;,U2i)B^iC2,(l>) (17) 

+ ^n:,n,2B^{n„<p)]. 

The poles in the Bl function are canceled as expected using the unitarity of the U, V and 
A^. The mass splitting is a function of Mi, M2, tan P, fi, M so ft and M. In the correct EW 
breaking parameter space, our numerical result for the mass splitting due to SUSY particles 
is that 

^^susY ~ O.Oig^ MeV (18) 

which is negligibly small compared to the pure EW corrections. 

B. The thermal relic density 

The thermal relic density fixes the WIMP mass. In the MDM, the relic species are 
{X^2,X^^,X°} and the coannihilation channels are: 

X^X^^AA,ff. (19) 

where A and / denote a EW gauge boson and the SM fermion, respectively. The mass 
splittings among them are very small compared to their masses. In the density calculation. 



such mass splittings are negligible. The relic particle thermal average cross section is 



fl: 



2 



7r«2 



{aAv){X^X^ ^ AA, ff)^^x 166 , (20) 

Matching to the relic abundance, the DM particle mass is determined to be M = 4.4 TeV 
without considering Sommerfeld corrections. 

Once SUSY is introduced in, a whole bunch of superpartners of those in MDM will 
present. SUSY breaking gives soft masses to the scalars and gauginos. Looking at X's, they 
have a mass M + Msoft- In general when M^oft < O.IM, eg. Msoft < 440 GeV for M = 4.4 
TeV, they have sizable effects on the relic abundance and must be included in the the relic 
species 6|, [Zl]. So now the coannihilation relic species are {(X^^,X^^,X*^), (X^^,X^^,X°)} 



and the coannihilation channels are: 

TX^ -^ AA, ff, GG, 

X!X^^ff,GA, (21) 

X^X^' ^ AA, ff, GG, 

where G denotes MSSM gauginos and / the superpatner of /. 

Furthermore, because M is also much larger than the electroweak scale, the physics of 
the above-mentioned coannihilation is basically supersymmetric and EW gauge symmetric, 
we can make unbroken EW symmetry and unbroken supersymmety approximation when 
calculating the thermal average cross section. Nevertheless , it is still a hard work. In 
terms of two-component fields, there are 24 gauge kinetic vertices and another 24 vertices 
involving superpatners. But actually we can have a useful and proper estimate for the 
relations between cross sections for the three kinds of processes in eq. (21). It is found that 
introducing in SUSY have nearly 4 times influence on the thermal average cross section. 
It will be shown in Fig. 2-8 a series of explicit examples by using two-component spinor 
techniques |8|] and the results reads: 

w(X+2x+2 -^ W3W3) = 16 X ^, (22) 

w(0*+V+' ^ W3W3) = 32 X ^, (23) 

va{X+^X+^ -^ CJ3CJ3) ~ 0, (24) 

w(0*+V^^ -^ ^3^3) ~ 0, (25) 



wcr(0+^0' ^ -)■ UJ3U3) ~ 16 X 



w(X+V^^ -> W^s^s) ^ 16 X 
w(X+20'"2 -^ Wscos) ~ 16 X 



37ra2 

2M2' 

37ra2 

8M2' 
37ra2 

8M2' 



(26) 
(27) 
(28) 



where v is the relative velocity in the lab frame, Eq. (24), (25) is the results of the p-wave 
suppression. 
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FIG. 2: The eight Feynman diagrams for X^'^X^'^ — > VF3VF3 
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FIG. 3: The three Feynman diagrams for 0*+20+2 _ > W3W3 
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FIG. 4: The four Feynman diagrams for X+2x+2 — )• cJstJs 
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FIG. 5: The two Feynman diagrams for 
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FIG. 6: The four Feynman diagrams for {^+^0 ^ 
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FIG. 7: The twelve Feynman diagrams for X 
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We have already calculate a complete series of thermal cross section, we end up with: 



fO", 



eff 



J2v(r^ 4t;a(X+2x+2 -^ W3W3 



(29) 



with considering the freedom for X+^, X+^ and 0^^,0 ~^ are all g = 2. This is very different 
with the dramatically influence from the higher representation such as SU{2)-5 in our prob- 
lem. Despite we make this conclusion from some specific examples, we think it is a common 
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result. 

Once the DM particle is as heavy as few TeV, the Sommerfeld effect due to hundreds GeV 
particles should be taken into consideration. In the MDM, the Sommerfeld enhancement 
effect due to SM particles, especially due to W and Z bosons, has been calculated, the 
factor is about 5 [9]. We expect approximately the same result in our case. While our case 
is A^ = 1 SUSY, all the SUSY particles should be included in the ladder diagram calculation 
in determining the potential which tells us the enhancement factor. However including the 
SUSY particles does not change the Sommerfeld effect essentially. In the extreme case of 
A^ = 4 SUSY, the extra symmetries just keep the gauge coupling from running. The potential 
itself has the same form as in the non-SUSY Yang-Mills case. In our A^ = 1 SUSY case with 
soft breaking, the logrithmic runing of the gauge coupling of the MDM is expected only 
mildly reduced, the Sommerfeld effect is then approximatelly the same. So it is reasonable 
to say that the Sommerfeld enhancement factor is about the same as that calculated in the 
MDM, M ~ VS X 2 X 4.4 ~ 19.7 TeV. 

C. Direct and Indirect signatures 

As for the direct DM detection rate, this SMDM is of the same order as the MDM. 
The DM particle interacts with quarks via loops. Although more heavy particles appear in 
the loops, the total cross section has the same order, o"si oc lO^^^cm^, which is within the 
sensitive of the current experiments, such as Super-CDMS and Xenon 1-ton lOl. Ill|. 

The DM annihilation in the galaxy may have observable signitures. The estimation of 
the cross section is also like that in the MDM 



12 



13|, the result has no order change. Note 



that SUSY does not change the Sommerfeld effect essentially, the predominant annihilation 
channel is still into EW W bosons, {av)ww ~ 10~^^cm^s~^. Because the DM mass M ^ 19.7 
TeV, we expect this model predict: (1) continuous rise e+/e+ + e~ spectrum up to about 20 
TeV; (2) fiat e+ + e~ spectrum up to M; (3) p/p flux has excess above the energy probed 



by PAMELA. SMDM has the solutions to the PAMELA and Fermi-LAT anomaly [ij, ll5|, 
the predictions in higher energies need further experimental data to verify. 
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IV. SUMMARY AND DISCUSSION 

We have made SUSY extension to the MDM model by introducing a complex scalar 
quintuplet as the superpartner of the fermion quintuplet. The neutral component of the 
fermionic 5-plet is still the DM particle as in the MDM model. Mass splittings among the 
fermionic 5-plets due to SUSY have been calculated in detail, they are found to be small. By 
considering new relic species and new coannihilation channels into the MSSM final states, 
the DM mass is estimated to be 19.7 TeV. 

The direct and indirect signals are basically the same as those in the MDM. Numerically 
the DM elastic scattering cross section with a nucleus is about 10^"^^ cm^ s^^ and the cross 
section of the predominant annihilation channel into W bosons is about 10~^^ cm^ s~^. 
SMDM predicts e+,e+ + e^ ,p spectrum in agreement with the previous PAMELA, Fermi- 
LAT data, p fiux has excess above the energy probed by PAMELA which need further 
experimental test. 

In the near future, suppose SUSY is discovered, say at LHC, it will be still a question 
if MSSM itself provides a DM particle, because R-parity as a discrete symmetry is still an 
assumption which is not as solid as gauge invariance and SUSY. It is plausible that R-parity 
is violated. In that case, it is still simple and interesting to have the DM via introducing 
SU(2)i high dimensional representations. 
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Appendix: the su(2)-n representation 

The SU(2)-n representation U^ ( n=2j+l ) is self-conjugate , when j is integer (n-odd), 
U^ is real , when j is half-integer (n-even), U^ is self-conjugate also but not real : 



xwx- 



w 



u^ = u^ 



j = 0, 1, 2, ... , X is symmetric 



U^ is not real j 



1 3 

2' 2' ■■• ' 



X is antisymmetric 



it is important in proving some identities including the SU(2) generators. 
The generators of the su(2)-n representation is: 
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(A.2) 



(A.3) 



(A.4) 



and have 



'Ta.'T^ -- rpZrp 
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(A.5) 



(A.6) 
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